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Abstract. We develop a statistical theory describing quantum- mechanical scattering of a particle by 
a cavity when the geometry is such that the classical dynamics is chaotic. This picture is relevant to a 
variety of physical systems, ranging from atomic nuclei to mesoscopic systems to microwave cavities; the 
main application here is to electronic transport through ballistic microstructures. The theory describes 
the regime in which there are two distinct time scales, associated with a prompt and an equilibrated 
response, and is cast in terms of the matrix of scattering amplitudes S. The prompt response is related 
to the energy average of S which, through the notion of ergodicity, is expressed as the average over an 
ensemble of similar systems. We use an information-theoretic approach: the ensemble of S-matriccs 
is determined by (1) general physical features— symmetry, causality, and ergodicity, (2) the specific 
energy average of 5, and (3) the notion of minimum information in the ensemble. This ensemble, 
known as Poisson's kernel, is meant to describe those situations in which any other information is 
irrelevant. Thus, one constructs the one-energy statistical distribution of S using only information 
expressible in terms of 5" itself without ever invoking the underlying Hamiltonian. This formulation 
has a remarkable predictive power: from the distribution of S we derive properties of the quantum 
conductance of cavities, including its average, its fluctuations, and its full distribution in certain cases, 
both in the absence and presence prompt response. We obtain good agreement with the results of the 
numerical solution of the Schrodinger equation for cavities in which the assumptions of the theory hold: 
ones in which either prompt response is absent or there are two widely separated time scales. Good 
agreement with experimental data is obtained once temperature smearing and dephasing effects are 
taken into account. 
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1. Introduction 

Scattering of waves by complex systems has captured the interest of physicists for a long time. For 
instance, the problem of multiple scattering of waves has been of great importance in optics 
Interest in this problem has been revived recently, both for electromagnetic waves || and for electrons 
Q, in relation with the phenomenon of localization, which gives rise to a great many fascinating 
effects. 

Nuclear physics, with typical dimensions of a few fm (lfm = 10 -15 m), offers excellent examples 
of quantum-mechanical scattering by "complex" many-body systems dating as far back as the 1930's 
when compound- nucleus resonances were discovered (5) . The treatment in these cases is often frankly 
statistical because the details of the many-body problem are intractably complicated. 

In the last two decades, electron transport in disordered metals has been intensively investigated 
|Q,|[0J|,^,|lO| , as has transmission of electromagnetic waves through disordered media HJlfl]. The 
typical size scale is 1/im for electronic systems and 1/im to 0.1m in the electromagnetic case. Because 
these are also examples of scattering in complex environments, where the character of the disorder is 
not exactly known, a statistical approach which treats an ensemble of disordered potentials is natural. 

Amazingly, features similar to those of these complex nuclear and disordered problems have also 
been found in certain "simple" systems. While the geometry of these systems is apparently very 
simple- quantum-mechanical scattering of just one particle by three circular disks in a plane, for 
instance- the classical dynamics is fully chaotic. Such systems have been studied by the "quantum 
chaos" community, in which the main question is how the nature of the classical dynamics influences 
the quantum properties |ll],[l^]. In contrast to the complex cases, these systems are amenable to 
exact (numerical) calculation; when the results are analyzed statistically, however, the results are 
closely related to those of the complex systems • Experimentally, two types of simple scattering 

systems have been studied in particular: electron transport through microstructures called "ballistic 
quantum dots", whose dimensions are of the order of lfim, and microwave scattering from metallic 
cavities, with typical dimensions of 0.1m. 

The "universal" statistical properties of wave-interference phenomena observed in systems whose 
dimensions span about 14 orders of magnitude turn out to depend on very general physical principles 
and constitute the central topic of this article. Although our main interest throughout the paper is 
electronic transport through ballistic chaotic cavities, or ballistic quantum dots, described in section 



1.2, we wish to emphasize the generality of the ideas involved, presenting first in section 1.1 their 
application in the field of nuclear physics- where some of them were first introduced- with a brief 
reference to microwave cavities at the end of that subsection. 



1.1. The atomic nucleus and microwave cavities 

One of the most successful models in nuclear physics, called the optical model of the nucleus, was 
invented by Feshbach, Porter and Weisskopf in the 1950's jl3],[l4|]. That model, which works very 
nicely over a wide range of energies, describes the scattering of a nucleon by an atomic nucleus- a 
complicated many-body problem- in terms of two distinct time scales: 

1. A prompt response arising from direct processes, in which the incident nucleon feels a mean field 
produced by the other nucleons. This response is described mathematically in terms of the average 
of the actual scattering amplitudes over an energy interval: these averaged amplitudes, also known as 
optical amplitudes, show a much slower energy variation than the original ones. 

2. A delayed, or equilibrated, response, corresponding to the formation and decay of the compound 
nucleus. It is described by the difference between the exact and the optical scattering amplitudes: 
it varies appreciably with energy and is studied with statistical concepts using techniques known as 
random matrix theory |l5| , |I(| . 

Just as in the field of statistical mechanics time averages are very difficult to construct and hence 
are replaced by ensemble averages using the notion of ergodicity, in the present context too one finds 
it advantageous to study energy averages in terms of ensemble averages through an ergodic property 
0,[8|Jl9@. 



The optical model not only works well in nuclear physics, but has also been applied successfully in 
the description of a number of chemical reactions, thus bringing us from the nuclear to the molecular 
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size scale 21 1. 



The connection of the above problems with the theory of waveguides and cavities was proposed 
very clearly by Ericson and Mayer-Kuckuk more than thirty years ago [0 : "Nuclear-reaction theory 
is equivalent to the theory of waveguides... We will concentrate on processes in which the incident 
wave goes through a highly complicated motion in the nucleus... We will picture the nucleus as a 
closed cavity, with reflecting but highly irregular walls." 

In fact, recent experiments with microwave cavities have shown features similar to those that had 
been observed in the nuclear case. Importantly, the "irregular walls" anticipated from the nuclear case 
are not necessary in order to see these features: the analogy between nuclear reaction theory and the 
theory of waveguides holds for simple smooth cavities as long as the corresponding classical dynamics 
is chaotic |f23| , p4] , p5f . This has been the focus of several recent experiments involving microwave 
scattering from metallic cavities |^6|,^7],^8 29 1 . These studies in simple chaotic systems were predated 



by extensive studies on scattering of microwaves by a disordered dielectric medium; particularly 
important are the experiments of A. Genack and collaborators |3(J. It is also interesting to note 
that the use of statistical concepts to analyze electromagnetic scattering by waveguides started quite 
independently of the connection to the complex scattering of nuclear physics, in the context of radio 
wave propagation pl| . 

1.2. Ballistic Mesoscopic Cavities 

The term mesoscopic system refers to microstructures in which the phase of the single-electron wave 
function- in an independent-electron approximation- remains coherent across the system of interest 
|Q,|U^J|,^,|lO| : this means that the phase-coherence length 1$ associated with processes that can change 
the environment- the other electrons, or the phonon field- to an orthogonal state exceeds the system 
size ||. This is realized in the laboratory with systems whose spatial dimensions are on the order of 
1/Ltm or less and at temperatures <1 K. In so-called ballistic cavities, or quantum dots, the electron 
motion is in addition practically ballistic, except for specular reflection from the walls: thus the 
elastic mean free path l e \ also exceeds the system size. In the most favorable material system, GaAs 
heterostructures, this condition can also be realized for cavities of size < lpm ||f7]j8|!|. 

Experimentally, an electrical current is established through the leads that connect the cavity to 
the outside and the potential difference across the cavity is measured, from which the conductance 
G is extracted. In an independent-electron picture one thus aims at understanding the quantum- 
mechanical single-electron scattering by the cavity, while the leads play the role of waveguides. It is 
the multiple scattering of the waves reflected by the various portions of the cavity that gives rise to 
interference effects. Three important experimental probes of the interference effects are an external 
magnetic field B, the Fermi energy cf, and the shape of the cavity: when these are varied the relative 
phase of the various partial waves changes and so the interference pattern changes. The changing 
interference pattern in turn causes the conductance to change; this sensitivity of G to small changes 
in parameters through quantum interference is called conductance fluctuations. 

The connection between scattering by simple chaotic cavities |23| , p4| , p5f and mesoscopic systems 
was first made theoretically |32| |. Subsequently, cavities in the shape of a stadium, for which 
the single-electron classical dynamics would be chaotic, were first reported in Ref. (3^]. More 
recently, several other types of structures have been investigated including experimental ensembles 
of shapes H|m|,||,|^,||,|9|,|^^,||,|3|,||] . Averages of the conductance, its fluctuations, and its full 
distribution were obtained over such an ensemble. 

It is the aim of the present paper to provide a theoretical framework to describe this physical 
situation. To this end we will set up a scheme similar to that explained in section 1.1 in connection 
with the scattering problem in nuclear physics. A complementary treatment based on semiclassical 
ideas has also been developed but will not be covered here |52|,fl5| . 

The paper is organized as follows. We start by presenting the general ideas of quantum-mechanical 
scattering by a cavity, introducing the scattering or S matrix of the problem (Sec. |2|), and then turn 
to treating an ensemble of systems in terms of an ensemble of S matrices (Sec. ^). Specific analytical 
results for the conductance are then presented, first in the absence of direct processes (Sec. |4|), and then 
with direct processes present (Sec. |^). We then compare our theoretical results with the numerical 
solution of the Schrodinger equation (Sec. [(]), and finally in Sec. compare with the experimental data 




Figure 1. The 2D cavity studied in the text. The cavity is connected to the outsi de v ia, L waveguides. 
The arrows inside the waveguides indicate incoming or outgoing waves as in Eq. (2.1). In waveguide / 
there can be Ni such incoming or outgoing waves: this is indicated in the figure by the amplitudes a l n , 
b l n , respectively, where n = l, ■ ■ -, JV; . 



that were already mentioned above. In the latter comparisons, a number of discrepancies are found; 
to reconcile theory and experiment, we realize the necessity of introducing the effect of processes that 
destroy the coherence of the wave function in the sample. Finally, the conclusions of our work are 
presented in Sec. g[ Some of the main results in this paper have appeared in condensed form in our 
previous publications jig^Q . 



2. The scattering problem 

2.1. Scattering waves: Definition of the S -matrix 

Consider a system of noninteracting electrons. Since we shall be dealing with cases in which spin-orbit 
coupling is negligible, we disregard the spin degree of freedom and just consider "spinless electrons" 
in what follows. We are interested in studying the scattering of an electron at the Fermi energy 
eF = H 2 kp/2m by the 2D microstructure shown schematically in Fig. [p. The microstructure consists 
of a cavity, connected to the outside by L leads, ideally of infinite length, that play the role of 
waveguides. The l-th lead (1 = 1,- ■ -,L) has width Wi. We are interested in the (scattering) solutions 
of the Schrodinger equation inside such a structure, with the ideal boundary condition that the walls 
of the cavity and leads are completely impenetrable: hence the wave function must vanish there. 

In each lead I we introduce a system of coordinates xi, yi, as indicated in Fig. [I]. The xi axis runs 
along the lead and points outwards from the cavity. The yi axis runs in the transverse direction and is 
tangential to the cavity wall and its continuation across the lead; yi takes on the values and Wi on 
the two walls of the lead. In lead I and for xi > we have the elementary solutions to the Schrodinger 
equation 

e ±tk " x 'Xn(yi) (2.1) 

where the positive (negative) sign is for outgoing (incoming) waves. Here the functions Xn(yi) are 
the solution of the transverse Hamiltonian which in the presence of a magnetic field may depend on 
kit The solution of the scattering problem consists in expressing the amplitude of the outgoing 

waves in terms of the incoming ones. 

In the absence of a magnetic field, the explicit problem can be simply stated; we now present this 
case in detail, noting that it can be generalized to the B^O case @. For B = 0, the functions Xn(yi) 



G 



arc 



(2.2) 



where is the "transverse wave number" . The functions Xn(yi) vanish on the two walls of the lead 
and form a complete orthonormal set of functions for the variable yi; i.e. 



(Xn\Xr, 



(2.3) 



This "transverse quantization" is a consequence of the boundary condition on the walls of the leads. 
Each possibility defined by the integer n is named a mode, or channel. The "longitudinal" wave 
number fe„ satisfies the relation 

[k£] 2 + [KPf = k* F . (2.4) 



If K$ < kp, then [fc^] 2 > 0, k$ is real, and the e ±lk ™ )xi occurring in Eq. (2.1) represent running 
waves along the leads: we thus have running modes or open channels. On the other hand, when 
K$ > kp, then [A;^] 2 < and kn is pure imaginary, thus giving rise to exponentially decaying 
waves along the leads: these modes are called evanescent modes or closed channels. If 



Ni < k F Wi/Tr < N t + 1, 



(2.5) 



there are Ni open channels in lead I. Very far away along the leads, i.e. as xi — > oo, only the 
contribution of the open channels contributes to the wave function. The most general form of the 
asymptotic wave function in lead / is thus the linear combination 



Ni 

E 

n=l 



e '^n X l 



(Hk^/m) 1 / 2 " (hk^/m) 1 / 2 



■b<P- 



(0 



XniVl)- 



Note that the normalization of the plane waves is such that they give rise to unit flux. 
We define the ^-dimensional vector 



i {l) ) T 

1 N l ) 



that contains the JVj incoming amplitudes in lead / (1=1, 
together, we form the vector 



,Ni). Putting all the (1 = 1 



(2.6) 



(2.7) 

..,£) 
(2.8) 



associated with the incoming waves in all the channels for all the leads. We can make similar definitions 
for the outgoing-wave amplitudes. The scattering matrix, or S matrix, is then defined by the relation 

b = 5a, (2.9) 

connecting the incoming to the outgoing amplitudes. In terms of individual leads we can write 



5' 



''li 
tai 



T21 



th\ tL2 



tlL 
t 2 L 

TLL 



(2.10) 



Here, ru is an Ni x Ni matrix, containing the reflection amplitudes from the iVj channels of lead I 
back to the same lead; ti m is an Ni x iV m matrix, containing the transmission amplitudes from the 
N m channels of lead m to the Ni channels of lead I. The S matrix is thus a square matrix with 
dimensionality n given by 



L 

E 

i=i 



(2.11) 
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Whil e we have carried out the analysis explicitly for B = 0, the main results, Eqs. (2.1), (2.6), (2.9), 
and ( 2.10 ), also hold in the presence of a magnetic field. 

Having chosen the unit-flux normalization for the plane waves in Eq. (|2.6| ), flux conservation 
implies unitarity of the S matrix jl 50 5l|j52]] ; i.e. 

S&=I. (2.12) 

In the absence of other symmetries, we have this requirement only. This is the unitary case, also 
designated in the literature as (3 — 2. In the presence of time-reversal invariance (TRI) (as in 
the absence of a magnetic field) and no spin, the S matrix, besides being unitary, is symmetric 
3l5lj5l,|5l: 



S = S J 



(2.13) 



This is the orthogonal case, also designated as (5=1. The symplectic case ((3 = 4), arising in the 
presence of spin and time-reversal invariance, will not be touched upon in this presentation. 



2.2. 



The conductance 



For a two-lead problem, which will occur most frequently in our analysis, the S matrix has the 
structure 



S 





tl2 ' 




' r t' ' 


. f 21 


T22 




t r> 



(2.14) 



In the particular case Ni=N2 = N, i.e. when the two leads have the same number of channels N, the 
four blocks r, t, r', t' are N x N and the S matrix is 2N x 2N. 

As we mentioned in the Introduction, we are interested in the electronic conductance of the 
microstructure. If we assume that the latter is placed between two reservoirs (at different chemical 
potentials) shaped as expanding horns with negligible reflection back to the microstructure, then the 
Landauer-Buttiker formula |]5ftp&p9| expresses the conductance G in terms of the scattering properties 
of the microstructure itself as 

G = jg, g = 2T, (2.15) 

where the factor 2 arises from the two spin degrees of freedom and the "spinless conductance" T is 
the transmission coefficient |i a f,| 2 summed over initial and final channels: 

T = Tr(tt t ). (2.16) 



2.3. Polar representation of the S matrix 

In the two-lead case with N\ — N2 = N, one can parametrize the S matrix in the so called "polar 
representation" as H,B3,bft|6[ 



5' 



Vi 










v 2 





v/7 





v 3 













VRW 



(2.17) 



Here, r stands for the ^-dimensional diagonal matrix of eigenvalues r a (a— 1, • • •, N) of the Hermitian 
matrix W . The Vi (i = l, • • -,4) are arbitrary N x N unitary matrices for (3 = 2, with the restrictions 

V3=V?, Vi = V% (2.18) 



for (3 = 1. It is readily verified that any matrix of the form (2.17) satisfies the appropriate requirements 
of an S matrix for (3 = 1,2. The converse statement, as well as the uniqueness of the polar 
representation, can also be proved using an argument similar to that in Ref. [55|]. 

In the polar representation, we can write the total transmission T of Eq. ( 2.1 6| ) as 



T 



(2.19) 



Thus the polar representation is natural for the study of conductance since it separates the magnitude 
of the transmission- the transmission eigenvalues {T a }~ from the irrelevant phase effects- the unitary 
matrices {vi}. 
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3. Ensembles of S matrices: An information-theoretic approach 

As we mentioned in the Introduction, we are interested in ensembles of systems, which will be 
represented as ensembles of S matrices endowed with a probability measure. For that purpose it 
is important to introduce first the notion of invariant measure for our S matrices: this we develop in 
the first subsection. The second subsection is devoted to the derivation of the probability measure for 
our ensemble of S matrices, using an information-theoretic point of view. 



3.1. The invariant measure 

The invariant measure is the measure which equally weights all matrices which satisfy the unitarity 
and symmetry constraints. Intuitively, it corresponds to the most random distribution consistent with 
the constraints, or in other words the one with the least information. Mathematically, such a measure 
is defined by requiring that it remain invariant under an automorphism of a given symmetry class of 
matrices into itself p3| , ^ |: 

d^ ] {S) = d^\S'). (3.1) 

For /3 = 1, the transformed matrix S' is related to S by 

S' = U SUo, (3.2) 



Uq being an arbitrary, but fixed, unitary matrix. Clearly, Eq. (3.2) is an automorphism of the set of 
unitary symmetric matrices into itself. For /3 = 2, 

S' = UoSVo (3.3) 

where Uq and Vq are arbitrary fixed unitary matrices. For (3 = 2, the resulting measure is the well 
known Haar's measure of the unitary group and its uniqueness is well known ]60|j6l|] . Uniqueness for 



/3 = 1 was shown in Ref. |53j. Using the invariant measures, Eqs. ( p.l[ ) - (3.3), as the probability 
measures for ensembles of S'-matrices defines the Circular Orthogonal and Unitary Ensembles (COE, 
CUE) for /3 = 1,2, respectively. 

Several explicit representations of the invariant measure are known, the classic one being in terms 
of the eigenphases and eigenvectors of the ^-matrix. For our purposes, the polar representation of Eq. 



(2.17) is of particular interest because of its connection to the conductance properties of the cavity. 
We thus consider the two-equal-lead case, N± = N2 = N, and express the invariant measure explicitly 
in this parametrization. 

We first recall a well known result from differential geometry. Consider the expression for the arc 
element 

ds 2 = }^g flv (x)6x fl 5x I/ , (3.4) 

written in terms of independent variables and the metric tensor g^ v (x). Assuming that ds 2 remains 
invariant under the transformation x ll =x^{x' 1 ,x l 2 , • ■ •), one can prove that the volume element 

dV = |detg(x)| 1/2 JJdx^ (3.5) 

A" 

remains invariant under the same transformation f62]| . 

We now go back to our random-5-matrix problem. We define the differential arc element as 

ds 2 = Tr[dSUS}. (3.6) 



This expression remains invariant under the transformations ( |3.2| ) and ( |3.3| ). Substituting for S the 
form ( |2.17 ), one can extract the metric tensor; applying Eq. (|3~5|), one then finds the invariant 



measure. This is done in [A|; the result is 

V« (S) = P<« ({t}) I] dr a J] dn(v® ) , (3.7) 
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where P^\{t}) denotes the joint probability density of the {r}, dfi(v^) is the invariant or Haar's 
measure on the unitary group U(N) EL and Cp is a normalization constant. From ^ we have 

P (1 \{r}) = Ct J] | r„ - r b | I] -1=, P (2) ({r}) = C 2 JJ | r„ - r 6 | 2 . (3.8) 

a<b c V c a<h 



Eqs. (3.7) and (3.8) explicitly specify the invariant measure. The factor involving the product over 
pairs gives the repulsion of the eigenvalues; notice that the repulsion is linear for the orthogonal case 
while quadratic in the unitary case, as typically occurs. 

3.2. The information-theoretic model 

3.2.1. The one- channel case To begin our discussion, consider first a physical problem that can 
be described by a 1 x 1 S matrix. This is the case, for instance, for a particle scattered by a ID 
potential that is nonzero in the region —a < x < 0, to which an impenetrable wall is added at x = —a: 
the particle then lives in the semiinfinite domain — a < x < oo. Another example is that of a 2D 
cavity, connected to the outside by only one lead that in turn supports only one open channel. From 
unitarity, 5* must be a complex number of unit modulus at every energy; i.e. S(E) = e l9i - E \ In the 
Argand diagram Re(S)-Im(S l ), S(E) is represented, for a given energy, by a point on the unitarity 
circle: that point is defined by the angle 9(E). As the energy changes, so does the representative 
point: this is what we may call, pictorially, the "motion" of S(E) as a function of energy; it resembles 
the motion in phase space of the representative point of a classical system as a function of time. 

We ask the following question: as we move along the energy axes, what fraction of the time do we 
find 9 lying inside a given interval d9? Let us call dP(9)=p(9)d9 that fraction. 

To answer this question, we first analyze how to construct energy averages. By this we always 
mean a local energy average, performed inside an interval / that contains many resonances and yet 
is small compared to an energy interval over which "gross structure" quantities showing a secular 
variation, like the average spacing A of resonances, vary appreciably. We then expect the dependence 
of the average in question on the center of the interval E$, its width /, as well as the particular 
weighting function used to define the average, to be weak. 

We devise an idealized situation: the argument E in S(E) is extended all the way from — oo to 
+oo, in such a way that local averages are everywhere the same as inside I: this idealization, which 
we call stationarity, will only represent well what goes on locally inside I in the actual system. We 
indicate an energy average of a quantity by placing a bar over it while an ensemble average is denoted 
by angular brackets. 

Ref. [|l9| chooses, as the weighting function, a Lorentzian. Using the fact that S(E) is analytic in 
the upper half of the complex energy plane (causality) , it shows that 

S* = [S]*, (3.9) 

i.e. the average of the k-th power of S coincides with the k-th power of the average of S. Thus, the 
quantity S, referred to as the optical S matrix in Sec. 1.1, plays a special role, in that the average 
of any power of S can be expressed in terms of it. Ref. [19|| then finds the answer to the above-posed 
question: the fraction of time p(9) spent by 9 in a unit interval around 9 in its journey along the 
energy axis is uniquely given by the expression 

-«-5^ S " W - (3 - 10) 

and depends only upon the average, or optical, S matrix S. This expression is also known as Poisson's 
kernel J54| . The conclusion is remarkable: it tells us that the system specific details are irrelevant, 
except for the optical S matrix. 

As an example, consider a ID (5-potential centered at x = and a perfectly reflecting wall at x= —a. 
An energy stretch containing 100 resonances starting from ka = 10, 000 (so that secular variations of 
gross-structure quantities can be neglected) was sampled to find the fraction of time that 9 fal ls in a 
certain small interval around 9 The result is compared in Fig. || with Poisson's kernel ( 3.10 ), 
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Figure 2. In the 5-potential model a stretch of energy containing 100 resonances, starting from 
fca = 10, 000, was sampled to find the fraction of time that 8 falls inajrnit interval around 9: the result 
is indicated by diamonds. The curve is a plot of Poisson's kernel (3.10), with the value of S extracted 
fron i the numerical data, in order to have a parameter-free fit. The agreement is excellent. (From Ref. 

where the value of the optical S was extracted from the numerical data themselves, so as to have a 
parameter-free fit. We observe that the agreement is excellent. 

Consider now a collection of systems, described by an ensemble of S matrices endowed with a 
probability measure. As an idealization, suppose we further construct S(E) as a stationary random 
function of energy, for — oo < E < oo. Then we know the conditions under which ergodicity, 
understood as equality of energy and ensemble averages (except for a set of zero measure), holds 
Q. Let us then assume that our ensemble is ergodic. 

The condition (3.9) arising from analyticity, together with ergodicity, implies the relation 



(s k ) = (sy 



(3.11) 



between ensemble averages, often called the analyticity- ergodicity (AE) requirement. The ensemble 
measure is thus uniquely given by 



dP, 



(S) 



(S)= P{S) (6)d6, p (s) (0) 



1 1 - 1 (S) I s 
2tt\S-{S)\< 



(3.12) 



once (S) is specified. The ensemble depends parametrically upon the single complex number (S), any 
other information being irrelevant. 



We note in passing that Eq. (3.11) implies that a function f(S) that is analytic in its argument, 
and can thus be expanded in a power series in S, must fulfill the reproducing property p4| , |66| 

(3.13) 



f((S)) = J f(S)dP {s) (S). 

It is because the probability measure appears as the kernel of this integral equation that it is called 
Poisson's kernel. 



3.2.2. The multi-channel case We now consider S matrices of dimensionality n, that can describe, 
in general, a multi-lead problem with n channels altogether, as explained in Sec. |^. The Argand 
diagram discussed above for n = 1 has to be generalized to include the axes ReSii, ImSii, Re5i2, 
11X1512, ■■■jReS'nn, ImS nn ] S is restricted to move on the surface determined by unitarity (SS^ = I) 
and, for j3=l, symmetry (S — S T ). 

We assume E is far from thresholds and recall that again S(E) is analytic in the upper half of 
the complex-energy plane. The study of the statistical properties of S is again simplified by idealizing 
S(E), for real E, as a stationary random-matrix function satisfying the condition of ergodicity. The 



same argument as in the ID case above shows that the AE requirement (3.11) is generalized to 



((S aibl ) ni ■ ■ ■ (Sa kbk T k ) - (S aibi r ■ ■ ■ (Sa kbk ) nk ■ (3.14) 
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Notice that this expression involves only S, and not S* matrix elements. Similarly, if /(S) is a function 
that can be expanded as a serie s of nonnegative powers of Sn, ■ ■ •, S nn (analytic in S), we must have 
the reproducing property (|3. 13 ) . 



Our starting point is the invariant measure dfip{S) that was introduced in the las t subsection. 



The average of S evaluated with that measure vanishes (shown explicitly in section IA ) , so that the 
prompt, or direct, com ponents described in the Introduction vani sh. I t is easy to check that the 
AE requirements ( 3.14| ) or, equivalently, the reproducing property (3.13), is satisfied exactly for the 



invariant measure. Ensembles that contain more information than the invariant one are constructed 
by multiplying the latter by appropriate functions of S. We relate the probability density pjm (S) to 
the differential probability through 

dP$(S)=pW(S)dw(S) (3.15) 

and require the fulfillment of the AE conditions. It was shown in Ref. |36| that, for n > 1, the AE 
conditions and reality of the answer are not enough to determine the probability distribution uniquely. 
However, it was shown that the probability density (Vp is a normalization factor) 



„ (q) v -l [d0t(J-(S) (<7)t)](/»n+2-/3)/2 

P(S){o) = Vg 7 , (3.1b) 

|det(/-S(S)T)|' 3 "+ 2 -' 3 



known again as Poisson's kernel, not only satisfies the AE requirements ( 1.14 ) |5j| , but the information 
X associated with it 



X\p] = j p {s) (S)ln P{s) (S)dfi(S) (3.17) 

is less than or equal to that of any other probability density satisfying the AE requirements for the 
(S) p6| ]. Notice that, for n = l, Eq. ( |3.16 ) reduces to (3.12). Thus the information entering 



same 



Poisson's kernel specifies : 1) General properties: i) flux conservation (giving rise to unitarity of the 
S matrix), ii) causality and the related analytical properties of S(E) in the complex-energy plane, 
and Hi) the presence or absence of time-reversal (and spin-rotation symmetry when spin is taken into 
account), that determines the universality class: orthogonal, unitary (or symplectic) 2) A specific 
property: the ensemble average (S) (=S under ergodicity), which controls the presence of prompt, or 
direct processes in the scattering problem. System-specific details other than the optical S are assumed 
to be irrelevant. 

The fact that for n > 1 AE and reality do not fix the ensemble uniquely is not surprising. In 
general (the lxl case being exceptional) we expect, physically, the matrix (S) to be insufficient to 
characterize the full distribution when, in addition to the prompt and equilibrated components, there 
are other contributions associated with different time scales [Q . Out of all possibilities, though, the 
information-theoretic argument selects the one where the prompt and equilibrated components and 
the associated optical S are the only physically relevant quantities. 

In addition to the completely general derivation of Poisson's kernel above, we present a concrete 
construction of this distribution following Ref. (6^j6^]. For the equilibrated part of the response, 
suppose there is an S-matrix 5*0 which is distributed according to the circular ensembles. For the 
prompt response, imagine a scattering process Si occuring prior to the response So- The total 
scattering is the composition of these two parts. Specifically, imagine bunching the L leads of the 
cavity into a "superlead" containing n incoming and n outgoing waves. Along the superlead, between 
the cavity and infinity, we connect a scatterer of the appropriate symmetry class described by Si. 
Since there are n incoming and n outgoing waves on either side of the scatterer, Si is 2n-dimensional 
and can be written 



Si = 



n t\ 



(3.18) 

The composition of the two scattering processes yields a total S 

S = n+f 1 (l-S r[)- 1 S h. (3.19) 
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One can prove |34|j6^j67| , |69| the following statement: the distribution of S is Poisson's measure 



( 3.16 ) with (S) — r\ if and only if the distribution of Sq is the invariant measure. That is, Eq. 



( 3.19 ) transforms between the problem with direct processes and the one without (for the one-energy 
distribution considered here). Also, one can show [54||39| that the distribution is independent of the 
choice of t\ and t\ , as long as they belong to a unitary matrix S\ . 

Note that throughout this work we use arguments which refer only to physical information 
expressible entirely in terms of the S'-matrix. An alternate point of view is to express everything in 
terms of an underlying Hamiltonian which one then analyzes using statistical or information-theoretic 
assumptions. These two points of view give, in fact, the same results: one can prove [[70H67|j6^1 that, 
for (S) = 0, a Gaussian Ensemble for the underlying Hamiltonian gives a Circular Ensemble for the 



resulting S. The argument was extended to (S) ^ in Refs. |67|j6^] using the transformation ( p,19[ ) 
above. 

4. Absence of direct processes 



In this case the optical matrix (S) vanishes and Poisson's kernel (3.16) reduces to the invariant 
measure: 

dP$ =0 = dfj,W(S). (4.1) 

We now derive the implications of this distribution for T, starting with the average and variance of T 
and then turning to its distribution. 

Averages of products of S: Weak-localization and conductance fluctuations 

Averages over the invariant measure of products of S matrix elements (invariant integration) can be 
evaluated using solely the properties of the measure, without performing any integration explicitly 
|^,^lj. We first discuss the unitary case, because it is simpler, and then the orthogonal one. 

4.. 1.1. The case [3 = 2 To illustrate the procedure, consider first the average over the invariant 
measure of a single S matrix element, to be denoted as 

(S aa ) ( o ] = I S aa d^\S). (4.2) 



Even though it is trivial to recognize that this average vanishes- since the invariant measure gives the 
same weight to S aa and to its negative- we present a more formal argument, which will be generalized 
later to more^ complicated averages. If U is an arbitrary but fixed unitary matrix, we define the 
transformed S as 

S = U°S. (4.3) 



Introducing (4.3) in ( f4.2| ) we have 

(S a 4 2) =X>°a' / S, a d^(5)^^(Mi 2) (4-4) 



where we have used the defining property (3.1) of the invariant measure and the definition of (S' a ' Q ) . 
In particular, if we take, as the arbitrary fixed matrix U 

[7° a ,=e^<W, (4.5) 

we find 

(S aa ) = e lBa (S aa ) ■ (4.6) 
Since this expression should hold for arbitrary 9 a7 we conclude that 

(S aa ) = 0. (4.7) 
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The above argument can be generalized to prove that 

(2) 



[<Sbl/3i ' ' 1 SbpPp] \S aiai ■ ■ ■ S aq a q 



o, 



(4.8) 



unless p — q and unless {ai, . . . , a p } and {&i, . . . , b p } constitute the same set of indices except for order, 
with the same condition for the sets {a±, . . . , a p }, {/3i, . . . , /3 p }. In particular, consider p = q— 1. Using 
the same argument as above, we find 



{Sb0S aa )l ' - £ U bb> { U aa') ( S b'(3S a , a ) [ Q ' , W 



(2) 



(4.9) 



For U given as in Eq. (|4.5| ), we have 



Sb0S* aa ) o = e^- e °^SbpS: a/ii 



(4.10) 



wh ich v anishes unless b = a. Had we defined S through right multiplication instead of left multiplication 
in (4.3), we would have concluded (3 — a. The only nonvanishing possibility in Eq. (4.9) is thus 



S a 



/ j | *-* aa' 



\s ala \ z ) , vc/°. 





(4.11) 



For instance, the choice of the matrix that produces a permutation of the indices 1 and 2 as U° yields: 



U° = 



1 

1 
1 







Sla\ 2 ) = (\S 2a \ 2 



(4.12) 



Similarly 

(|5 lQ | 2 ) o = --- = (|5„ Q | 
The final result for the average intensity follows from unitarity 

(2) 1 



a=l 



1 



So 



(4.13) 



(4.14) 



As an application, we calculate the average conductance when we have a cavity connected to the 
outside by means of two leads supporting Ni and N 2 open channels (so that n = N\ + ./V2): 



Nt N 2 
a=l b=l 



(2) 



mm 



Ni + N 2 



1 1 

~N~! + iV2 



(4.15) 



which is the series addition of the two conductances N\ and N 2 . (The superscript on the angular 
brackets indicates (3 — 2.) 

With similar arguments one finds [fn|: 



S 12 \ Z \S; 



34 



5'l2| 2 I ^13 1 2 



(2) 


(2) 



ri 2 - 1 



1 



(2) 



n(n + 1) ' 
2 

n(n + 1) ' 



(4.16) 
(4.17) 
(4.18) 
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Here, 1,2,3,4 stand for any quartet of different indices, so that n in each case must be large enough to 
accomadate as many indices as necessary. 

As an application, we calculate the second moment of the conductance as 

Ni N 2 

,2\(2) 



a,c 

of the condu' 

[var(T)]< 2 > = 



2 u ,2 
cd\ 



(2) 



N?Nl 



ii 



Q ,c = i M= i (N 1 + N 2 y-i 

The variance of the conductance is then given by 



(jvi + w 2 ) 2 [(iVi + N 2 y - 1 

Notice that in the limit Nx, N 2 -> oo with Ni/N 2 = K fixed 

A^ 2 



[var(T)]( 2) 



(4.19) 



(4.20) 



(4.21) 



a constant which depends only on the ratio Ni/N 2 and on no other details of the cavity. Since this 
limit corresponds to increasing the width of the waveguides and so of the full system, the fact that the 
result is constant is the analog of the so-called universal conductance fluctuations (UCF) well-known 
for quasi-lD disordered systems ||,[|. I n particular, for K=l, var(T) — > 1/16, slightly less than the 
quasi-lD value of 1/15. 

4-1-2. The case (3=1 Just as above, we first illustrate the procedure through the average over the 



invariant measure of a single S matrix element [Eq. (4.2) written for (3=1] which vanishes trivially. 



We introduce the transformed S, again a unitary symmetric matrix, through 

S = £/ S(£/°) T (4.22) 

Substituting (4.22) in the integral definition of (S a b)o and using the defining property of the invariant 
measure (3.1) we have 



In particular, if we take as the arbitrary fixed matrix U° the one given by Eq. (4.5), we find 



(S ab ) Q = e^+^> (S ab ) Q =► (5 ab ) = 

since 9 a , 9b are arbitrary. 

Just as for (3 = 2, the above argument can be generalized to prove that 



S aibl ■ ■ ■ S ap b p [SaiPi ' ' ' ^a q Pq\ 



(1) 



= 0, 



(4.23) 



(4.24) 



(4.25) 



unless p = q and unless {ai, b±, . . . , a p , b p } and {ai, (3i, . . . , a p , (3 P } constitute the same set of indices 
except for order. In particular, for p = q=l we find 



S<tbS a p) = E U aa' U bb> (Uaa'Upp') (S a >b> S a , p,\ , VJ7°. 
a'b' a' 3' 



For instance, for U given as in Eq. (4J5), we have 



(4.26) 



(4.27) 



which vanishes unless {a,6} = {a,/3} or {/3,a}. As a particular case, take a = b = a = (3=1 and, as U°, 
a matrix with the structure 

U° 2 ••■ 



u° 



u& C/ 2 °2 o 
1 





(4.28) 
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We find 



s "l 2 ) = E U ^u° w {ul,ul v ) (s a , v sl, v 

a'b' 

+ ^la'^lfe' (^Sy^la') \Sa f b'Sb'a' 



a'b' 
— 

= [|^i|V|^ 2 r](|5 11 | 2 ) o + 4|[/ 1 1 | 2 |C/ 1 ° 2 | 2 (|ft 2 | 2 ) o . (4.29) 
Squaring the unitarity relation jt/^] 2 + |t^i 2 | 2 = l and substituting in Eq. (4.29) we finally obtain 

(|^ 11 | 2 )^-2(|5 12 | 2 ^ 1) . (4.30) 

This result is very important. It states that time-reversal invariance (TRI) has the consequence that 
the average of the absolute value squared of a diagonal S-matrix element is twice as large as that of 
an off-diagonal one, under the invariant measure. By unitarity, the specific value of each one of these 
averages is given by 

(!) 2 /, l2 \« 1 



S a ^b\ 



(4.31) 



Just as in the above case (3 = 2, we calculate as an application the average conductance when our 
cavity is connected to the outside by two leads with N\ and N 2 open channels (n = N\ + N 2 ): 

Ni N 2 



(i) 



mN 2 



a=l b=l 



N!+N 2 + l 



(4.32) 



Here, the extra 1 in the denominator as compared with Eq. (4.15) is the weak-localization correction 
(WLC), a symmetry effect resulting from TRI. We can rewrite Eq. ( 4.32 ) separating out the WLC 
term as 

(1) NiN 2 N ± N 2 
(T)o = JVi + N 2 ~ {Ni + N 2 ) {Ni +N 2 + 1)' (433) 

In particular, for Ni= N 2 = N and for TV — > oo, corresponding to a large system, the WLC term tends 
to the universal number —1/4. 

Using a similar procedure, one finds the results {n being again the dimensionality of the S matrix) 



S 12 | 2 |S 34 | 2 



Si 2 \ 2 \S 13 \ 2 



W_ n + 2 

o n{n + \){n + 3) 

(i) 1 



Sis 



(1) 



o n{n + 3) ' 
2 



(4.34) 
(4.35) 
(4.36) 



o n{n + 3) 

A comment like that made right after Eq. (4.18) applies here as well. Just as in Eq. ( 4.19| ), we now 



find for the second moment of the conductance 

. , i , _ N X N 2 [NiN% (,Vi + N 2 + 2) + 2] 



<^>o 



(Ni + N 2 ){Ni +N 2 + l)(iVi + N 2 + 3) 
and for its variance 

[var(T)]«- 2^(^ + 1)^ + 1) 



K 2 



{N 1 + N 2 ) {Ni +N 2 + i f {^ + N 2 + 3) {K + l) 4 



(4.37) 



(4.38) 



where in the limit N±, N 2 — > 00 with N\/N 2 — K, a fixed number. Note that in this universal limit, 
the variance here is exactly twice as large as for (5 = 2, Eq. (4.21), another result of time-reversal 
invariance. For the particular case K=l, the limiting value of the variance is 1/8. 
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4-2. The distribution of the conductance in the two- equal-lead case 

In the last section we characterized the conductance of a cavity through the first two moments of 
its transmission. If the distribution of the conductance is Gaussian, this is sufficient to characterize 
the full distribution. In fact, it can be shown that in the large-size universal limit, N — > oo, the 
distribution of the conductance is indeed Gaussian |72| . 

In general, however, the probability density of T will not be Gaussian, and it is of interest, then, 
to derive results for this density. For this purpose, the polar representation of Sec. |2.3| is particularly 
useful since the conductance is directly related to the {r} whose joint probability distribution we 



know. Specifically, the distribution of the transmission T of Eq. (2.19) can be obtained by direct 
integration of the P^({t}) of Eq. (O): 



} ( T ) = J S(T-Y,Ta)P m ({Ta})l[dT a . (4.39) 

a a 

We consider a few examples below. 

4-2.1. The case N = 1 In this case we have only one r , that we may call r, so that T — t, and 



< T < 1. Eq. then gives 



W W(T) = ^=, wW(T) = l. (4.40) 

For f3—l, we thus have a higher probability to find small T's than T ~ 1: this is clearly a symmetry 
effect, a result of TRI that favors backscattering and hence low conductances. 

4.2.2. The case N = 2 Now T = t 1 + r 2 , and < T < 2. In || we show that 

|T, < T < 1 

w (1) (T)={ ' (4.41) 



and 



{T-2^/T^\) ) 1<T<2 



wW(T) = 2[1- |1-T|] 3 (4.42) 



For /3= 1, notice the square-root cusp at T =1. We find, once again, a higher probability for T < 1 
than for T > 1 to occur. On the other hand, for /3 = 2, w(T) is again symmetric around T=l. 

4.2.3. The case N = 3 Now T=n + r 2 + t 3 , and < T < 3. For (3 = 2 one finds 

w {2) (T)={ , 7KS ,-, 88rr 101o7 , 2 ^ 100 ^ 3 (4-43) 



7 -T— 1818T 2 + 1836T 3 
-1035T 4 + 324T 5 - 54T 6 + fT 7 - |T 8 , 1< T< § 

and the distribution is symmetric about T=1.5. As mentioned above, w^'(T) gradually approaches 
a Gaussian distribution. 

4-2.4- Arbitrary N In this case it is straightforward to obtain the dependence of the tail of the 
distribution in the region < T < 1. In this region the constraint that r < 1 does not enter; a 
calculation given in [b| shows that 

w^(T)ocT^ n2 / 2 -\ (4.44) 
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5. Presence of direct processes 



In order to treat cases involving direct processes, as in billiards in which short paths produce a prompt 
response, we now need Poisson's kernel ( 3.16| ) in its full generality. We discuss below some analytical 
results for the distribution of the spinless conductance T in the case of a cavity connected to the 
outside by means of two leads supporting one open channel each (N\ = N2 = 1), giving rise to a 
2-dimensional S matrix. There is only one r in this case and it is its distribution that we seek, since 
T = t. While the expressions that we derive are somewhat cumbersome, they are used for comparing 
to numerical results in section |^ where plots of several examples are displayed. 
We write the optical S-matrix S, a subunitary matrix, as 



x w 
z y 



(5.1) 



where the entries are, in general, complex numbers. 
5.1. The case (3 — 2 

From Eq. ( 3.1 6| ) we write the differential probability for the S matrix as 



<)(s)= [det(/-Ss¥ ^ )(s)| 
\det(I-SS J )\ 2n 



where 



V 

We are interested here in the case n = 2. 



d^(S) = l. 



In the polar representation ( 2.17 ) with n = 2, the S matrix has the form 



S = 





e if} 



e l ~< 
e' 



the invariant measure ( p.3| ) being [see Eqs. (3.7), (3.8)] 

dad(3dr/dS 



dfia(S) = dr- 



(2nf 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



1) As a particular case, suppose the optical S matrix is diagonal, so that there is only direct 
reflection and no direct transmission: in Eq. 



Substituting Eqs. flEUD, (|5j), flS^ in (|J)we find 



5.1) we choose w = z = 0. 



dP^(T,<p,fl>) 



(l-X 2 ) 2 (l-Y 2 



drdipdip 



|(e-^ + X % /T — 7) (e-^ ^ YVT~^ r) - XYt\ (2tt) z 
where ip = a + 7, tp = (3 + 5, X = \x\, Y= \y\. The distribution of the conductance T is thus 



(2) 
,J XY 



(T) = {l-X 2 ) 2 (l-Y 2 



(5.6) 



(5.7) 



where 



where 



(e-^ + X^/l~^T) (e-^ - Y^/l~T) - XYT\ 



) . denotes an average over the variables ip, -0. The result is (0 < T < 1) 

(2) fT , _ K A - B(l — T) + C(l - T) 2 + D(l - Tf 
[E-2F(1-T) + G(1-T) 2 } 5/2 

K = (l- X 2 ) 2 {1 - Y 2 ) 2 , A = (1 - X 4 Y 4 )(l - X 2 Y 2 ) 



(5.8) 
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B = (X 2 + y 2 )(l - 6X 2 Y 2 + X 4 Y 4 ) - 
C = (1 + X 2 Y 2 )(6X 2 Y 2 - X 4 - Y 4 ) - 
D = (X 2 + Y 2 )(X 2 -Y 2 ) 2 , E = 

f = (i + x 2 y 2 )(x 2 + r 2 ) - 4x 2 y 2 , 



4X 2 Y 2 (1 + X 2 Y 2 ) 
4X 2 Y 2 (X 2 + Y 2 ) 
[1 - X 2 Y 2 ) 2 

G= (X 2 -Y 2 ) 2 . 



This result reduces to 1 when X =Y = 0, as it should. A particularly interesting case is that of 
"equivalent channels", i.e. X = Y, in which the above expression reduces to 

(2) (T) n {l-X 4 ) 2 + 2X 2 (l + X 4 )T + AX 4 T 2 

w x,x\ T ) = ( 1_ x > — 



{l-X 4 ) 2 

[(1 - X 2 ) 2 + AX 2 T] 5/2 



(5.9) 



The structure of this result is clear if we notice that 



(2) 

,J x,x 



(0) 



l + X 2 



l-X 2 



> 1, 



(2) 

,J xx 



(1) 



(l-X 2 ) (l+X 4 
(l + X 2 ) 3 



< 1, 



(2) (2) 

and hence w^yo) > w x x (l), so that small conductances are emphasized, as expected, because of 
the presence of direct reflection and no direct transmission. 

2) The case of only direct transmission and no direct refle ctio n is obtained by setting x = y = in Eq. 
(5T). The conductance distribution is obtained from Eq. ( |5.8j) with the replacement X — > W=|to|, 
Y— ► Z =\z\, T — ► 1- T. In the equivalent-channel case we now obtain a conductance distribution 
that emphasizes large conductances. 

3) The case of a general optical S matrix, Eq. (5.1), has also been worked out and the result 
expressed in terms of quadratures: because of its complexity, it will not be quoted here. 



5.2. The case /3 = 1 

This case is more complicated than that for f} = 2 and we have only succeeded in treating analytically 
some particular cases. Take, for instance, a diagonal S, i.e. w — z = in Eq. ( |5.l| ). With the same 
notation as above, we find 

W ^ Y (T) = \(l-X 2 )(l-Y 2 )\ 3/2 ^= (5.10) 



(e- i{ f + Xy/T^T) (e-^ - FyT^T) - XYT 



3 / 



a result that has to be integrated numerically. When X = Y — 0, the distribution (5.11) reduces to 
1/2-v/r, as it should. It is interesting to notice that for Y = the above result can be integrated 
analytically, to give 

(l - x 2 ) 3 ^ 2 

w%(T) = V 7 2 f\ (3/2,3/2;l;X 2 (l-T)) , (5.11) 

2-Fi being a hypergeometric function [Q . 



6. Comparison with numerical calculations 

The information-theoretic approach that we have been discussing is expected to be valid for cavities in 
which the classical dynamics is completely chaotic, a property that refers to the long time behavior of 
the system. It is in such structures that the long time response is ergodic and equilibrated, and so one 
can expect that maximum entropy considerations will play a role. In this section we examine particular 
cavities numerically in order to determine to what extent the information-theoretic approach really 
holds. The structures that we consider are all "billiards"- they consist of hard walls surrounding a 
cavity with constant potential- with two leads. We start by considering particularly simple structures, 
then treat structures in which the absence of direct processes is assured, then move on to structures 
having particularly obvious direct processes. 
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Figure 3. The magnitude of (a) the quantum correction to the classical conductance and (b) the 
conductance fluctuations as a function of the number of modes in the lead N. The two asymmetric 
stadium-like billiards shown at the right were used; the average of the results for the two cavities is 
shown. The numerical results for B = (squares with statistical error bars) and B ^ (triangles) are 
compared to the predictions of the COE (points, dotted line) and CUE (points, dashed line). The 
agreement is poor. For both billiards, 25 energies were sampled for each value of iV, and for non-zero 
field BA/<fco = 2, 4 where A is the area of the cavity. 



6.1. Simple structures 

In the "quantum chaos" literature- the study of how quantum properties depend on the nature of 
the classical dynamics in a system- several billiards are used as standard examples of closed chaotic 
systems. The two most studied are the Sinai billiard- the region enclosed between a square and a 
circle centered in the square- and the stadium billiard- two half-circles joined by straight edges. The 
classical dynamics in these two billiards is known to be completely chaotic. For a test case open 
system, then, it is natural to take one of these billiards and attach two leads. The open stadium 
billiard was studied previously for this reason ||^,|3^,|7^,[76],(77| . Here we directly compare results for 
this system to the predictions of the information-theoretic approach. 

The numerical methods used in these calculations are covered in detail in Ref. |74[ . Briefly, the 
procedure consists of the following three steps. First, discretize the Hamiltonian onto a square mesh 
using the simplest finite-difference scheme. Solving the Schrodinger equation is then reduced to solving 
a set of linear equations. Second, find the Green function at the desired energy from one lead to the 
other using a recursive procedure and outgoing-wave boundary conditions. This procedure essentially 
uses the sparseness of the finite-difference matrix to efficiently solve the linear equations. Third, note 
that the transmission amplitude can be obtained from this Green function by simply projecting onto 
the transverse wave-functions in the lead. The main parameter in these calculations, ka, is the size 
of the mesh compared to the wavelength. In the results shown here, ka is always less than 0.8 and 
ka < 0.5 in most cases. For these values, the anisotropy of the Fermi surface is small; the non- 
parabolicity of the dispersion is larger, but does not concern us here since we treat transport at a 
fixed energy. 

Two simple quantities to calculate both numerically and theoretically are the average and variance 
of the conductance. The analytic results for the invariant measure are in section 4.1. Fig. || shows 
the numerical results for the two asymmetric open stadium shown on the side- asymmetric half- 
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Figure 4. Statistics of the eigenphases of the S-matrix for the second stadium shown in Fig. g 
with N = 9. (a) Cumulative distribution of the eigenphase density (solid line) compared to the CE 
(dashed), (b) Cumulative distribution of the difference between nearest-neighbor phases for both B = 
and BA/ifio = 4 compared to the COE (dotted) and CUE (dashed). The spacing is normalized to 
the mean separation, (c) Variance of number of phases in an interval L for both B = (squares) and 
BA/ 4>q = & (triangles) compared to the COE (dotted) and CUE (dashed). All three statistics agree with 
the prediction of the circular ensembles- constant density, a spacing distribution given by the Wigner 
surmise, and a logarithmically increasing variance- despite the poor agreement for the transmission in 
Fig. |. 



stadiums are used in order to avoid the complications of reflection symmetry. The top panel shows 
the deviation of the average transmission from the classical value of the transmission. This classical 
value was obtained numerically by tracing trajectories through the cavity. For fully equilibrated 
scattering, the classical value is N/2 where N is the number of channels in each lead, but in Fig. || 
^classical 1S (0-58) for the upper (lower) cavity. The bottom panel shows the variance of T. 
While the numerical results are similar in magnitude to the predictions, the agreement is clearly not 
very good. 

Before proceeding with our discussion of the conductance in these cavities, we step back to perform 
the most common test of random matrix theory. In the context of closed systems, it is natural to 
consider the statistics of the energy levels. The degree to which these statistics agree with the Wigner- 
Dyson statistics derived from random matrix theory is often used as the prime indication of the validity 
of the theory for a given system. In the context of the S-matrix, the analog is to look at the statistics 
of the eigenphases: in the large N limit their statistics is also Wigner-Dyson Jll],[l6|]. In fact, studies 
of eigenphases of chaotic scattering systems were carried out prior to any interest in the conductance 
[ 23 1 24 1 ■ The statistics of the eigenphases can be characterized by three representative quantities 
[15,16|. First, the mean density of eigenphases measures the uniformity of the system; for the circular 
ensembles (CE) it is constant. Second, the nearest neighbor distribution highlights the repulsion at 
short scales; it is approximately the Wigner surmise for the CE. Third, the variance of the number of 
phases within a certain range L, denoted S 2 (L), indicates the rigidity of the spectrum at large scales; 
it grows logarithmically with L for the CE. 

These three quantities are shown for the simple open stadium in Fig. ^. The agreement with 
the predictions of the CE is good for all three, in contrast to the results for the conductance above, 
despite N not being very large. Since the conductance involves the transmission coefficient which is a 
property of the wavefunctions, this indicates that the distribution of the eigenvectors is more sensitive 
to deviations from the CE than the distribution of eigenphases. Thus the eigenphase statistics cannot 
be taken as a definitive indication of the validity of the CE: it is perfectly possible to have excellent 
eigenvalue statistics while having poor eigenvector statistics JT8| . 

Returning to the properties of the conductance, we believe that the deviations from the CE seen 
in the numerics (Fig. ||) are caused by the presence of short paths in these simple structures. This 
means that the response is not fully equilibrated. The two most obvious types of short paths in these 
structures are the direct paths between the leads and the whispering gallery paths- those that hit only 
the half-circle. Short paths will be included in the analysis in section 6.3 below. One way to minimize 
the effect of these short paths is to make the openings to the leads very small so that the probability 
of being trapped for a long time increases. Presumably the CE will apply to any completely chaotic 
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Figure 5. The magnitude of (a) the quantum correction to the classical conductance and (b) the 
conductance fluctuations as a function of the number of modes in the lead N. The numerical results for 
B = (squares with statistical error bars) agree with the prediction of the COE (points, dotted line), 
while those for B^O (triangles) agree with the CUE (points, dashed line). The six cavities shown on 
the right were used; the average of the numerical results is plotted. Note that each cavity has stoppers 
to block both the direct and whispering gallery trajectories. For non-zero field, £?j4/</)o = 2, 4 where A 
is the area of the cavity. (After Ref. [M.) 

billiard in the limit that the openings are very small (the number of modes in each lead should remain 
constant). However, such a structure is difficult to treat numerically, except in the very small N limit, 
which, in fact, we discuss in section |6.3| below. 



6.2. Absence of direct processes 

In order to compare with the predictions of the circular ensembles, we wish, therefore, to study 
structures in which the most obvious direct paths are absent. To this end, we have introduced 
"stoppers" into the stadium to block both the direct and the whispering gallery trajectories: examples 
are shown in Fig. In order to study the statistical properties, the conductance as a function of 
energy is calculated. Because the energy variation is on the scale of ?i7 esc (the escape rate from the 
cavity) |^ , ^2| , it is much more rapid than the spacing between the modes in the leads (fivp/W). 
Thus many independent samplings of the conductance at a fixed number of modes may be obtained. 
In addition, we vary slightly the position of the stoppers so as to change the interference effects and 
collect better statistics. The numerical results in Fig. || used 50 energies for each N (all chosen 
away from the threshold for the modes) and the 6 different stopper configurations shown; the classical 
transmission probabilities for these cavities ranged from 0.46 to 0.51 with a mean of 0.49. In addition, 
for non-zero magnetic field two values were used. We see that the agreement with the CE is now very 
good for both the mean and the variance, both for B = and for non-zero B. This supports our view 
that the deviations in the simple structure Fig. || are caused by short paths. 

While the agreement of the mean and variance of the conductance with the CE gives a strong 
indication of the validity of the information-theoretic model for real cavities, a much more dramatic 
prediction of the model is the strongly non-Gaussian distribution of T for a small number of modes. 
The analytic results for the CE were given in section |4.2| . These are compared to the numerical results 
for N = 1, 2 in Fig. ^|, using the same data as in Fig. ||. Note that the data is consistent with a 
square-root singularity in the case N = 1 B = and with cusps in the two N = 2 cases. Thus we see 
that even for this much more stringent test, the agreement between the behavior of real cavities and 
the CE is excellent. 
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Figure 6. The distribution of the transmission intensity at fixed N = 1 or 2 in both the absence and 
presence of a magnetic field, compared to the analytic COE and CUE results. The panels in the first 
row are histograms; those in the second row are cumulative distributions. Note both the strikingly non- 
Gaussian distributions and the good agreement between the numerical results and the CE in all cases. 
The cavities and energy sampling_Doints used are the same as those in Fig. ^ for By^O, BA/<f>o=2, 3, 
4, and 5 were used. (After Ref. |46[ .) 



6.3. Presence of direct processes 

We now want to look at more general structures than those used in the last section for comparing to 
the circular ensemble results. In particular, we will remove the stoppers that blocked short paths, and 
compare the numerical results with the predictions of Poisson's kernel, following Ref. |4§| . We will 
not, however, consider the most general structure: Poisson's kernel is expected to hold in situations 
where there are two widely separated time scales, a prompt response and an equilibrated response. 
Thus we will study structures where we expect this to be true. Since we have only obtained explicit 
results in the case N=l (see section ||), we will study this case numerically as well. 

We have computed the conductance for several billiards shown in Fig. (?]. Statistics were collected 
by (1) sampling in an energy window larger than the energy correlation length but smaller than the 
interval over which the prompt response changes, and (2) using several slightly different structures. 
Typically we used 200 energies in kW/n <E [1.6,1.8] (where W is the width of the lead) and 10 
structures found by changing the height or angle of the convex "stopper" . Note that the stopper here 
is used to increase statistics, not block short paths. As in the absence of direct processes, since we 
are mostly averaging over energy, we rely on ergodicity to compare the numerical distributions to the 
ensemble averages of the information-theoretic model. The optical S-matrix was extracted directly 
from the numerical data and used as (S) in Poisson's kernel; in this sense the theoretical curves shown 
below are parameter free. 

We first consider the billiard shown in Fig. at low magnetic field (BA/<fio = 2 where A is the area 
of the cavity; by low magnetic field we mean that the cyclotron radius r c is much larger than the size 
of the cavity (r c — 55 W shown to scale). In this case w(T) is nearly uniform [Fig. 0(a)], and (S) is 
small because direct trajectories are negligible in this large structure. We thus obtain good agreement 
with the invariant-measure prediction of a constant distribution ( 4.40) ). 



In order to increase (S) we make one of three changes: (1) introduce potential barriers at the 
openings of the leads into the cavity (dashed lines in structures of Fig. [?]), (2) increase the magnetic 
field, or (3) extend the leads into the cavity. In the first case, the barriers are chosen so that the bare 
transmission of each barrier is 1/2 (determined by calculation in an infinite lead). They cause direct 
reflection and skew the distribution towards small T [Fig. 0(b)]. Since the reflection from the barrier 
is immediate while the transmitted particles are trapped for a long time, one has two very different 
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Figure 7. The distribution of the transmission coefficient for N = l in a simple billiard (top row) and 
a billiard with leads extended into the cavity (bottom row) . The magnitude of the magnetic field and 
the presence or absence of a potential barrier at the entrance to the leads (marked by dotted lines in 
the sketches of the structures) are noted in each panel. Cyclotron orbits for both fields, drawn to scale, 
are shown on the left. The squares with statistical error bars are the numerical results; the lines are 
the predictions of the information-theoretic model, parametrized by an optical S-matrix extracted from 
the numerical data. The agreement is good in all cases. (From Ref. pj|.) 

response times. Second, the large magnetic field (BA/<po = 80) corresponds to r c just larger than the 
width of the lead (r c = 1.4 W). The field increases one component of the direct transmission- the 
one corresponding to skipping orbits along the lower edge- and skews the distribution towards large 
T [Fig. 0(c)]. Third, extending the leads into the cavity increases the direct transmission in both 
directions and also skews the distribution towards large T [Fig. @(d)]. We have done this rather than 
consider a smaller cavity since in our case the equilibrated component is trapped for a long time, 
yielding a clear separation of scales. 

In each of these cases, the numerical histogram is compared to the information-theoretic model 
(solid lines) in which the numerically obtained S is inserted. In panels (b)-(d) the curve plotted is the 
analytic expression of Eq. (5.S) and the corresponding one for direct transmission. Note the excellent 
agreement with the information-theoretic model. 

Since the long-time classical dynamics in each of the three structures (a), (b) and (d) is chaotic, 
these results show that a wide variety of behavior is possible for chaotic scattering, the invariant- 
measure description applying only when there is a single characteristic time scale. 

In case (c), the dynamics is not completely chaotic because of the small cyclotron radius, and 
so one would not expect the circular ensemble to apply. In Ref. jf6| we found that increasing the 
magnetic flux through the structure beyond a few flux quanta spoiled the agreement with the circular 
ensemble; we now know that a nonzero (S) is generated and that the present model describes the 
data very well. The excellent agreement found here with a flux as high as 80 suggests extending the 
analysis to the quantum Hall regime. 

By combining several of the modifications used above, different (S) and so different distributions 
can be produced. First, by using extended leads with barriers at their ends, one can cause both 
prompt transmission and reflection; the result is shown in Fig. 0(e). Finally, increasing the magnetic 
field in this structure produces a large average transmission and a large average reflection. The 
resulting w(T), Fig. 0(f), has a surprising two-peak structure: one peak near T — 1 caused by the 
large direct transmission and another near T—l/2. For cases (e) and (f), four intervals of 50 energies 
each were treated independently (since the four intervals show slightly different (S)'s ) and the four 
sets of data were then superimposed |frjj|| . Even in these two unusual cases, the prediction of the 
information-theoretic model is in excellent agreement with the numerical results. 
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7. Comparison with experimental data: Dephasing effects 

The random S'-matrix theory of quantum transport through ballistic chaotic cavities was seen in the 
last section to be in good agreement with numerical simulations for structures in which the assumptions 
of the theory are expected to hold: that agreement includes the average conductance, its variance and 
probability density. 

The ultimate test of the theory should, however, be comparisons with experiment. A naive 
comparison indicates poor agreement. For instance, for Ni — N2 = N = 2, both the theoretical 
weak-localization correction (WLC) and variance are larger than the experimental results |3S{] . The 
theoretical WLC (including spin factor) is —2/5 = — 0.40 from Eq. (|4.33 ), while the experimental one 



is 0.31. The theoretical var( 5 ) is 4/15 w 0.27 for /3 = 2 [Eq. (4.20)] and 72/175 w 0.41 for (3=1 



Eq. ( t.3S| )], giving the ratio [var(g)]p^ / [var(g)][, 2 ^ =54/35 « 1.54. The experimental results for the 



variance, on the other hand, are ~ 0.015 for B ^ and ~ 0.034 for B = 0, giving the ratio 2.27. In 
addition, the measured probability density |3q] is close to a Gaussian distribution, which differs from 



the prediction of Eqs. (4.41), (4.42). More recent experimental data for N= 1 show a distribution 



which is clearly asymmetric for B = 0, favoring low conductances as required by weak localization, but 



the asymmetry is by no means as pronounced as that of Eq. (4.40) which has a square-root singularity 
at the origin. 

To reconcile these discrepancies, we must realize that inherent in the discussion of the previous 
sections is the assumption that one can neglect processes which destroy the coherence of the wave 
function inside the sample and neglect energy smearing caused by non-zero temperature. Even if the 
phase-breaking length 1$ is larger than the geometrical size of the cavity, for sufficiently narrow leads 
the electron may spend enough time inside the sample to feel the effect of phase-breaking mechanisms. 
The discussion of these effects and their relevance for the description of the experimental data is the 
subject of the present section. 

We simulate the presence of phase-breaking events through a model invented by M. Buttiker [|S0|], 
where, in addition to the physical leads 1,2 attached to reservoirs at chemical potentials fi\, /12, a "fake 
lead" 3 connects the cavity to a phase-randomizing reservoir at chemical potential ^3. The idea is that 
any particle which exits the cavity through lead 3 is replaced by a particle from the reservoir; since 
the replacement particle comes from a reservoir, it is incoherent with respect to the exiting particle, 
hence the phase-breaking. Requiring the current in lead 3 to vanish determines the two-terminal 
dimcnsionless conductance is then found to be 

T23T31 



G/(e 2 /h) 



T 9 



T32 + T31 



(7.1) 



In this equation, the factor of 2 accounts for spin explicitly and TJy is the transmission coefficient for 
"spinless electrons" from lead j to lead i. This expression for g has a natural interpretation: the first 
term is the coherent transmission and the second term is the sequential transmission from 1 to 2 via 
3. In terms of the S matrix, is 

Nt Nj 2 

r « = E E|C bj | ■ ( 7 - 2 ) 

Qi = l b j = l 

where Ni is the number of channels in lead i; S^ hj is the matrix element in the position a,, bj of the 
ij block of the S matrix, rows and columns in this block being labeled from 1 to Ni and from 1 to 
Nj, respectively. The total number of channels will be designated by 

N T = J2 N i- (7-3) 

i 

Occasionally, we shall use the notation N$ to designate the number of phase-breaking channels N$ in 
the fake lead 3. These N^ channels are physically related to the phase-breaking scattering rate 7^ via 
the relation N^/^Ni + N 2 ) ~ "f<f>/jcsc, where 7 CSC is the escape rate from the cavity. This "fake lead" 
model has been used for various studies of the effect of phase-breaking in mesoscopic systems Q . 
We now make the assumption that the total Nt x Nt scattering matrix S obeys the distribution 



( [4.lD given by the invariant measure [|47|]82| 1 . Through this assumption, the effect of the third lead is felt 
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somehow uniformly in space rather than at any given point. This even-handed statistical treatment 
makes the "fake lead" approach more physically reasonable. A further generalization to the case of 
non-ideal coupling of the fake lead to the cavity is studied in Ref. ( |g3|) but will not be presented 
here. In the following, we confine the analytical discussion to the large N$ limit, then present the 
results of numerical random-matrix theory calculations, and finally compare again to experiment. 



7.1. Large 

We find below the average and variance of the conductance when 3> 1 p7| . As we shall see, we 
need, for that purpose, the average and the covariance of the transmission coefficients Ty introduced 
above. 



From Eqs. (4.14) and (4.31), the average of (i ^ j) is given by (as usual, /3 = 1, 2) 
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(7.4) 

(7.5) 
(7.6) 
(7.7) 



where the ~ sign refers to t he situation A3 ^> 1 while Ni, N 2 — 0(1). Turning to the covariance, using 
Eqs. (4.16)-(4.18) and (Q, one finds for /3 = 2 



N 2 y [N T S lk 5 3 i - N k N T 8 3 i - N l N T S lk + N k N t ] 



(7.6 



Likewise, using Eqs. (4.34)-(4.36) for (3 = 1 yields 



N T (N t + 1) (N T + 3) 
x {NiNj {N T + 1) (Nt + 2) (8 ik 8 fl + S u 8 jk ) + 2N i N k S ij 8 kl 
+2N i N k N l 8 ij + 2N i N J N k 5 kl + 2N t N T (N T + 1) S ijkl + 2N l N N k N l 
- (N T + 1) [2NiNi8 ijk + 2N i N k 8 ij i + 2N i N j (8 iM + 8 jkl ) 
+N i N j N l (8 lk + 8 jk ) + NiNjN k (8 U + 8 jt )}}. 



(7.9) 



Here, a 8 with two or more indices vanishes unless all its indices coincide, in which case its value is 1. 
These expressions are valid for arbitrary Ni, Nj, and also for i=j if Ta is interpreted as the reflection 
coefficient Ra from lead i ba ck t o itself. 

The conductance of Eq. ( |7.l| ) can be expressed as 
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where we have written 
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For convenience, we use interchangeably a bar or the bracket (• • -)^' to indicate an average over the 
invariant measure. 



From Eqs. (7.8) and (7.9) we find, for the variance of Ty (i ^ j) 
2 \(2) _ N Z N, (N T ~ N t ) (N T ~ N 3 ) 
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and 
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, \ (i) _ NiNj [(N T + 1 - Ni) (N T + 1 - Nj) - (N T + 1) + N.N,} 



(7.12) 
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N T (N t + 1) (N T + 3) 

For two leads only, i.e. for = 0, these last two expressions reduce to Eqs. ( 4.20| ) and ( 4.38 ), 
respectively. On the other hand, for N$ ^> 1 and i = l,2 



2\ W 



1/2 



T 



i?, 



OI Nl 



(7.14) 



Thus, ST23/T23 and ST13/T13 are small quantities in Eq. (7.10), and we can make the expansion 



g = 2(T 21 +ST 21 )+ 2 23 ^ 
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The average of this conductance for N<j, ^ 1 is obtained by using Eq. ( 7.14 ): 
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Substituting Eqs. (7.5)-(7.7), we have 
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which gives the weak-localization correction 
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Turning to the variance, we subtract the average conductance from Eq. (7.15) and obtain to lowest 
order in l/N^ 
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and obtain, for (3 = 1,2, to lowest order in l/N^ 
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(7.19) 
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Squaring this last expression and averaging, we obtain the variance of the conductan ce i n terms of 
the (8Tij8Tki)o- We must now substitute the variances and covariances given in Eqs. (7^8) and (7.9) 
to obtain our final result: 
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Figure 8. Magnitude of quantum transport effects as a function of the number of channels in the leads, 
Ni = N2 =N, for N$ = 0, 2, 4, and 8. (a) The weak-localization correction, (b) The variance for the 
orthogonal case (B = 0). (c) The variance for the unitary case (nonzero B). Open symbols are numerical 
results (20, 000 matrices used, statistical error is the symbol sizel. Solid lines are interpolation formulae. 
Solid symbols are experimental results of Refs. [ pi (squares), |35| (triangles), and Jj^ ] (circles) corrected 
for thermal averaging. The introduction of phase-breaking decreases the "universality" of the results 
but leads to good agreement with experiment. (From Ref. [EtJ.) 



The ratio of variances for f3=l and (3 — 2 is 



[var( 5 )]{, 2) 
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(7.22) 



We observe from this last equation, first, that this ratio is independent of to leading order, and, 
second, that the ratio of variances is larger than 2, and as high as 3 for N 2 = Ni = N = 1. For 
comparison, in the absence of phase-breaking and for N 2 — Ni, that ratio lies between 1 and 2. 



7.2. Arbitrary 

In order to evaluate effects of phase-breaking when is not large, we numerically evaluate the 
random-matrix theory, generating random Nt x Nt unitary or orthogonal matrices and computing g 
from Eq. ([7j]). Fig. || shows the weak- localization correction and the variance of the conductance as 
the number of modes in the leads is varied for several fixed N^. This result is relevant to experiments 
at fixed temperature in which the size of the opening to the cavity is varied. Though Sg and varg are 
nearly independent of N in the perfectly coherent limit- the "universality" discussed above- phase- 
breaking channels cause variation. Thus the universality can be seen only if <g; N; otherwise, the 
behavior is approximately linear. 

The results of Sec. ||, which correspond to — 0, and those of the last section, obtained for 
3> 1, suggest an approximate interpolation formula for arbitrary values of N$. For the case 
N 2 = Nx = N, Ref. 



Jf7| proposes 
N 



(7.23) 
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Figure 9. Probability density of the conductance in the orthogonal (first column) and unitary (second 
column) cases for N = I (first row) and N = 2 (second row) . Increasing the phase-breaking from zero 
(dashed lines, analytic) to Nj, = 1 (plus symbols, numerical) to N$ = 2 (squares, numerical) smooths 
the distribution. (From Ref. Et| . ) 
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for the variance. These interpolation formulae are compared with the numerical simulations in Fig. 
|§|: the agreement is good, the only significant deviation being for N = l and small N$. 

In addition to the mean and variance, the probability density of the conductance, w(g), can be 
evaluated. Fig. ^| shows the distribution in the weak phase-breaking regime with a small number of 
modes in the lead. From the variance and mean, we know that the phase-breaking will narrow the 
distribution and make it more symmetric because the weak-localization correction goes to zero. Not 
surprisingly, the phase-breaking in addition smooths the distribution, and the extreme non-Gaussian 
structure in w(g) is washed out. 



7.3. Experiment 

Now that we have a tool for evaluating the effects of phase-breaking, we return to a comparison of the 
information-theoretic model with experiments. First, for the mean and variance of the conductance, 
in Fig. H are also shown as black dots the results of three experiments from Refs. |}4 35 3^] . Since the 



theoretical results take into account through the effect of finite temperature insofar as dephasing 
is concerned, without including thermal smearing, the experimental variance has been corrected to 
compensate for this latter effect |Q. In the case of Ref. |38| measurements of all three quantum 
transport quantities were made, and so this data can be used to test the consistency of the theory. 
Notice that a value of = 4-8 allows the simultaneous fit of the WLC and the variance for both j3 = 1 
and 2. This indicates that the theory and experiment are indeed in good agreement. 

Second, the full distribution of the conductance has been measured in the N = 1 case and compared 
to the random-matrix theory including phase-breaking p4| . The results are shown in Fig. |l0|. In the 
B = case, the experimental distribution is clearly not Gaussian and is skewed towards small values 
of <?, as expected. On the other hand, the distribution is not nearly as dramatic as the completely 
coherent theory suggests. An important recent development in relation with the fictitious lead model 
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Figure 10. Experimental conductance distributions for both B = (open circles) and 40 mT (filled 
circles) for a 0.5 /im 2 device at 100 mK with N = l. These are compared to theoretical curves for both 
zero temperature (dashed lines) and non-zero temperature (solid lines). Though the effect of non-zero 
temperature is substantial, the measured distribution at B = is clearly not Gaussian. Inset: pattern 
of gates defining the dot. (After Ref. [Q.) 

described above has been carried out in Ref. ||S3|| , with the motivation of describing a spatially more 
uniform distribution of phase-breaking events: the fictitious lead is considered to support an infinite 
number of modes, each with vanishing transmission, allowing a continuous value for the dephasing 
rate. The experimental results are analyzed in Ref. |Q with this improved model, and the result 
is also shown in Fig. 
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It is seen that the shape of the conductance distribution is reproduced 
well with a dephasing rate obtained independently from the WLC. This comparison is strong support 
for the validity of the theory. The numerical value of the variance at various temperatures shows a 
discrepancy, though; in particular, the observed ratio of variances for (3 — 1 and P — 2 is considerably 
larger than that given by the model: this is as yet unexplained. 



8. Conclusion 



In this paper we developed a statistical theory aimed at the description of the quantum-mechanical 
scattering of a particle by a cavity, whose geometry is such that the classical dynamics of the 
system is chaotic. We studied, as our main application, the electronic transport through ballistic 
microstructures, in an independent-particle approximation. 

The theory, which was developed in the past within nuclear physics to describe the scattering of a 
nucleon by a nucleus, describes the regime in which there are two distinct time scales, associated with a 
prompt and an equilibrated response. The prompt response is described in terms of the energy average 
S— also known as the optical S matrix- of the actual S matrix. Through the notion of ergodicity, 
S is calculated as the average (S) over an ensemble of similar systems, represented mathematically 
by an ensemble of S matrices belonging to the universality class in question: orthogonal, unitary or 
symplectic. In addition, the ensemble satisfies the analyticity (causality)-ergodicity requirements and 
the optical (S) evaluated over it has a specified (matrix) value. The ensemble discussed in the text 
is the one that carries minimum information, or maximum entropy, satisfying these conditions. It is 
thus meant to describe those situations in which any other additional information is irrelevant. In this 
procedure one constructs the one-energy statistical distribution of S using only the above physical 
information- expressible in terms of S itself- without ever invoking any statistical assumption for the 
underlying Hamiltonian which never enters the analysis. 



30 



From the resulting S'-matrix distribution, known as Poisson's kernel, properties of the quantum 
conductance have been derived: its average, its fluctuations, and its full distribution in certain cases, 
both in the absence and presence of direct processes. We obtain good agreement with the results of 
the numerical solution of the Schrodinger equation for cavities in which the assumptions of the theory 
hold: either ones in which prompt response is absent (section ^), or ones in which there are two 
widely separated time scales (section ||) . As for the comparison with experimental data, agreement 
was found once temperature smearing and dephasing effects were taken into account, at least within 
a phenomenological model. 

The effect of time-reversal symmetry- determining the universality class (3- has been seen to be 
of fundamental importance. One other symmetry was not touched upon in this article, because of 
lack of space: spatial reflection symmetries. While such symmetries are not relevant for disordered 
systems- the traditional subject of mesoscopic physics- they are possible in ballistic systems in which 
one can control the scattering geometry. How such symmetries affect the interference contribution 
to transport was studied in Refs. |^,^{|. For reflection symmetries, S is block diagonal in a basis 
of definite parity with respect to that reflection, with a circular ensemble in each block (if (S) =0). 
The key point is that the conductance may couple the different parity-diagonal blocks of 5 1 , and thus 
the resulting quantum transport properties are a nontrivial generalization of the circular ensemble 
results. These effects are discussed in Ref. ||||6| for structures presenting a "left-right", "up-down", 
"four-fold" , and "inversion" spatial symmetry and compared with the results of numerical solutions of 
the Schrodinger equation. From an experimental point of view, it would be nice to have a confirmation 
of our results, both in microstructures as well as in microwave cavities, the latter ones being perhaps 
simpler to handle. 

We should point out that the theory presented here is not applicable when there are other relevant 
time scales in the problem. One example is the case of a disordered quasi-lD system, where the 
diffusion time across the system is also an important time scale. In the past, this extra piece of 
information was taken into account explicitly by expressing the transfer matrix for the full system as 
the product of the transfer matrices for the individual slices that compose the system |87| , |5^ , lS8| ,fi||7|]l0|l . 
An example of an intermediate situation is that of a finite number of cavities connected in series. 
Whether more than two time scales are physically relevant for single cavities as well is not clear at 
present. It seems lik ely that the deviations of the numerical results from the theory for the simple 
structures in Section 6A are for this reason J89|. Also, it is conceivable that additional time scales are 
at the root of the discrepancies discussed in Ref. |^5| in relation with the WLC. An extension of the 
maximum-entropy theory presented above to include other time scales is not known at present. 

We also stress that the theory presented here is applicable to the statistics of functions of the S 
matrix at a single value of the energy. The joint statistical distribution of the S matrix at two or more 
energies has escaped, so far, an analysis within the philosophy described above (some aspects of the 
two-point problem have been studied assuming an underlying Hamiltonian described by a Gaussian 
ensemble as in Refs. [ p0|j9l| ]). An approach coming close to that philosophy was initiated in Ref. 
p2[ for the simplest quantity of a two-point character: the statistical distribution of the time delay 
arising in the scattering process- a quantity involving the energy derivative of S- motivated by the 
study of the electrochemical capacitance of a mesoscopic system. It was found that another piece of 
information is needed for the description: the statistical distribution of the if -matrix of resonance 
widths. Ref. deals with lxl S'-matrices; subsequent work derived the distribution of time delays 



for N channels ]93j. In a further development, Ref. ]94| finds a transformation that relates the fc-point 
distribution of the n-dimensional 5-matrix for the case (S) ^ to that for (S) = 0, thus relating the 
problems in the presence and in the absence of direct processes. 

Finally, we remark that the information-theory described in the text makes use of the standard 
Boltzmann-Shannon definition of information and entropy. Other definitions have been presented in 
the literature: see, for instance, Ref. and the references contained therein. What would be the 
physics behind the use of other definitions of entropy and how our results would be modified is not 
known at present. 
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A. Appendix: Evaluation of the invariant measure 



We evaluate dS in the polar representation (2.17) and then the arc element ds 2 of Eq. ( |3.6| ) keeping V 
and W independent, as is the case for = 2. We use this same algebraic development and set W = V T 



in the proper place to analyze the /3 = 1 case. Differentiating S of Eq. (2.17) we obtain 



dS = (dV)RW + V(dR)W + VR(dW) = V [(5V)R + dR + R{SW)} (A.l) 
where we have defined the matrices 

5V = VUV, 8W = {dW)W\ (A.2) 

which are antihermitian, as can be seen by differentiating the unitarity relations V^V = I, WW^ =1. 
The arc element of Eq. (3.6) is thus 

ds 2 = Tr{ [R T (SV) ii + dR T + (SW^ R T ] ■ [{5V)R + dR + RSW}} 

= Tr[R T (SV)^(SV)R + R T (SV)UR + R T (SV)^R(SW) 

+ {dR T )(6V)R + {dR T )dR + (dR T )R{5W) 

+{5W)' i R T {6V)R + (6W) f R T dR + (SW) 1 R T R(SW)]. 

The various differentials occurring in the previous equation can be expressed as 



(A.3) 



dR=- 
2 



SV 



6W 



dr/^fp dr/y/? 
dr/y/r -dr/y/p 

v\dv\ 
v\dv2 



SV! 

Sv 2 



(dv 3 )vt 
{dvi)v\ 



Sv 3 
8v± 



where we have used the abbreviation 

p = 1 — T. 

We now calculate the various terms in (A.3). The first plus ninth terms give 

4 

Tr [(SV)1{SV) + {5W)\8W)] = Ti^(5v{)Uv,. 



(A.4) 
(A.5) 
(A.6) 

(A.7) 
(A.8) 



The second plus fourth terms give 

Tr [R T (8V)UR+(dR T )(SV)R] (A.9) 

= r Tr[ v /p(<5u 1 ) t <y/5+ yft dyf? + yp? {Svrf dyf? + y/p{Sv 2 )^ dy/p + h.c] = 0, 

where h.c. stands for Hcrmitian conjugate. We have used the identity 

Tr [D(5v^D' + D'(Sv)D] = 0, (A.10) 
where D and D' denote any two real diagonal matrices and 5v an antihermitian one. Using the identity 



Tr \(Sv l )' 1 D (Svj) d] = Tr \d (Svrf D {6v t ) 



(All) 
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we find for the third plus seventh terms 

Tr [R T (SVy R(SW) + (5W^R T (5V)R] 

= 2Tr[Vp {SvJ VP (^s) + ^ (Sv 2 )^ V7 (dv 3 ) 

+^(5v 1 )^(Sv i ) + ^(Sv 2 )^(6v 4 )]. (A.12) 
The fifth term gives 

Tr(^)(^) = i^^l!. (A.13) 

A a TaPa 

Finally, the sixth plus eighth terms give 
Tr [(dR T )R(SW) + h.c] 

= ^Tr [-(d,T)5v 3 + {dr)6v 3 + {dr)6v 4 - {dT)Sv 4 + h.c] = 0. (A.14) 
Substituting these expressions in (A. 3), we find 

ds 2 = 2Tr 1 1 ( 5v J + {dT lrp T) + VP (Sv 3 ) 

+ ^ (<^ 2 ) f (Svs) + ^ {Svi) 1 (Sv 4 ) + VP (SvJ V? (5v 4 ) } • (A.15) 

The antihcrmitian matrix Svi can be expressed as 

Svi = 5di + iSsi, (A. 16) 

where 5ai is real antisymmetric and Ssi is real symmetric. Substituting in the expression for ds 2 and 
rearranging terms, we find 

ds " = E ( E ((^)„„) 2 + 2pa [(S Sl ) aa (Ss 3 ) aa + (Ss 2 ) aa (Ss 4 )J 



+ 2r a l(5s 2 ) aa (S S3 ) aa + (5 Sl ) aa (5s 4 ) aa ] + g^M 

+ 2 E IE ((^)a h ) 2 + E ((^) ab ) 2 (A-17) 

a<fc L i=l i=l 

+ 2 VP^b[( fe i)a6 ( 5 «3) o6 + (<M ah (<M a6 + {5si) ab (Ss 3 ) ab + (5s 2 ) ab (Ss 4 ) ab ] 
+ 2^T b [(5a 2 ) ab (Sa 3 ) ab + (Jai) o6 (<5a 4 ) afc + (£s 2 ) ab (5s 3 ) a6 + {5 Sl ) ab {Ss 4 ) ab }} . 

A.l. The case (3 = 1 

In this case, v 3 = vf, v 4 = v 2 so Sv 3 = deltavi) T , 5v 4 = deltav 2 ) T , and hence 

<5a 3 = — 5ai, Ss 3 = 6s\, Sa 4 = —5a 2 , Ss 4 = 8s 2 . (A. 18) 

Substituting in (A. 17), we have 

ds 2 = 2j2{((Ssi) aa ) 2 + ((Ss 2 ) aa f + 2r a (5 Sl ) aa (5s 2 ) aa 

a 

{dT a f \ 



+ 



+ ■ 



+4E {((^l)a fc ) 2 + ((^)a fc ) 2 + ((^OaJ 2 + ((^) ab ) 2 
a<fc 

+2 v ^¥ b '[(Ss 1 ) ab (Ss 2 ) ab - (Jai) o6 (<5a 2 )J 

+ VP^ [((^i) ah ) 2 + ((M ab ) 2 - ((<foi) o6 ) 2 - ((^ 2 ) afc ) 2 ] } , (A.19) 
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-ds 2 = (1 + Pa) \((Ssi) aa ) 2 + ((Ss 2 ) aa ) 2 } + 2r a (5 Sl ) aa (Ss 2 ) aa + 



+2J2 {a + vwd [((<^u 2 + (o^ur 



(dTg) 2 

^TaPa 



(A.20) 



a<b 



+ (1 - y/Wb) [((^i) a6 ) 2 + (O^U) 2 ] + 2^T b [(S Sl ) ab (Ss 2 ) ab - (5 ai ) ab (Sa 2 ) ab ]j. 

In (A.20), (Ssi) . {Ss 2 ) aa and r a (a=l, • • -,N) contribute N independent variations each, while 
(^i)afc' (^2)067 {^ a i)ab^ a ^)ab ( a < &) contribute iV(iV - l)/2 each, giving a total of 



= 3N + 4 



N(N-1) 



2N 2 + N, 



(A.21) 



which is the correct number of independent parameters for a 2iV-dimensional unitary symmetric 
matrix. 

The metric tensor appearing in Eq. (3.4) has a simple block structure, consisting of 1 x 1 and 
2x2 blocks along the diagonal, as follows. There are TV 2-dimensional blocks with rows and columns 
labeled (5 Sl ) aa , {5s 2 ) aa : 



(Ssi) aa (5s 2 ) aa 

1 + Pa T a 
T a 1 + Pa 



(A.22) 



and N 1-dimensional blocks l/Ap a T a labeled by dr a , giving, altogether, the contribution 

N r-\ , \2 2 N 1 

n (i±£Lzi_ n I (A . 23) 

1 ^Pa 'a t ' a 

a—1 a—1 

to det g. There are N(N — l)/2 2-dimensional blocks with rows and columns labeled (6si) b , ($s 2 ) ab , 

(Ss%) ab {Ss 2 ) ab 

(A.24) 



( 5s l)ab 



1 + y/PaPb \jT a Tb 
\Jr a Tb 1 + y/PaPb 



■2. 



and also N(N — l)/2 2-dimensional blocks with rows and columns labeled {Sai) b , (8a 2 ) ab , 

( 6a l)ab ( Sa 2)ab 



( 6a l)ab 
( Sa 2)ab 



1 - y/PaPb -y/T a T b 
-\fT a Tb 1 - ^JPaPb 



(A.25) 



■2. 



giving, altogether, a contribution 



l + V(l-r«)(l-7i) 



1- V(l-Ta)(l-TJ,) 



2-r a -r b + 2V(l-r )(l-r & ) 2 - t q - r b - 2^(1 - r a )(l - r 6 ) 



= (r - r b ) z 



(A.26) 



Multiplying (|A.23[) . 
result given in Eq. (3.S 



and (A.26) and taking the square root, as required by Eq. ( |3.5| ), we find the 
)■ 
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A.2. The case (3 = 2 

We go back to Eq. (A. 17). We can write the single summation in that equation (except for its last 
term) as 

£{[(*«i)«„ + + i(^) aa + NJ' ~ 2r„ [(S Sl ) aa - {6s 2 )J [(Ss 3 ) aa - (8s,)J} 

a 

= £ [(fa a ) 2 + {5y a ) 2 - 2T a {5z a )(5x a - Sy a - 8z a )] , (A.27) 



where we have defined the combinations 
Sx a = (S Sl ) aa + (Ss 3 ) aa 

Sy a = (Ss 2 ) aa + (Ss 4 ) aa 



(A.28) 
(A.29) 

= (^Oaa ~ 0**0., . (A.30) 

Notice that in ds 2 the 4iV quantities {8si) aa , (8s2) aa , {Ss 3 ) aa , (Ss4) aa appear only through the 3N 
combinations Sx a , Sy a , 5z a : these quantities, together with the dr a contribute 47V independent 
variations. The (Ssi) ab and the (8ai) ab for i = 1, • • ■, 4 and a < b contribute 4 • N(N — l)/2 each, so 
that we have a total of 4N 2 variations, which is the correct number of independent parameters for a 
2/V-dimensional unitary matrix. 

In terms of independent variations we can thus write the ds 2 of Eq. (A. 17) as 

ds2 = E I [( fe «) 2 + (^) 2 - ^ a (Sz a )(6x a - Sy a - 5z a )) + 1 



(A.31) 



+2E E(W 2 +E(WJ 2 

o<h I z=l i=l 

+2VpIp6 [(^ai)afc (<M a6 + (<M ab (<5a4) ab + {Ssi) ab {Ss 3 ) ab + (Ss 2 ) ab (Ss 4 ) ab \ 
+2^ [(*02) o6 (<M a & + (<5ai) ab (<5a 4 ) ah + (Mat, (<5s 3 ) Qb + (S Sl ) ab (Ss 4 ) ab }}. (A.32) 



The metric tensor appearing in Eq. (3A) has a simple block structure, consisting of 1 x 1, 3 x 3 
and 4x4 blocks along the diagonal, as follows. There are TV 3-dimensional blocks with rows and 
columns labeled Sx a , Sy a , 5z a : 





Sx a 


Sy a 


Sz a 




Sx a 


1 





-T a 


(A.33) 







1 


Tq 




Sz a 


-To 


T a 


2r a 





and N 1-dimensional blocks l/2p a r a labeled by dr a , giving, altogether, the contribution 

N 



n 2r Q (l -Ta) = 
. 2t«(1 - r a ) 

a— 1 v 7 



(A.34) 



to det g. There are N(N — l)/2 4-dimensional blocks with rows and columns labeled (6si) b , (8s2) ab , 
( 5s 3) ab , (Ssi) ab (for a < b) 



( Ss l)ab ( Ss 2)ab ( 5s 3) ab ( Ss *)c 



(<fol)«6 
(^2) ab 
(<^) ab 
(^4) ab 




(A.35) 
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Figure Bl. Region of integration for finding the distribution w(T) in the N = 2 case. 

and also N(N — l)/2 4-dimensional blo cks w ith rows and columns labeled (Sai) ab , (Sa2) ab , (Sa3) ab , 
(<5a4) ab and identical to the matrices of (A. 35) giving, altogether, a contribution 

N 4 



(r a - n) 



(A.36) 



Multiplying ( A. 34 ) and ( A.36 ) and taking the square root, as required by Eq. (3.5), we find the 
result given in Eq. (3.S). 



B. Appendix: The distribution of the conductance in the absence of direct processes 



We derive here the two-channel conductance distribution of Eqs. ( 4.41 ) and ( 4.42 ) and the behavior 
of that distribution for arbitrary N in the range < T < 1, as given by Eq. ( 4.44 ). 

For the two-channel distribution we have to perform the integral (4.39) in the two-dimensional 
space ti, T2, inside the square < n < 1, < t 2 < 1, indicated in Fig. Bl. 



T = Ti + r 2 , 



n - t 2 



The change of variables 
(B.l) 



will be found advantageous. As shown in Fig. Bl, when < T < 1 the variable r varies in the interval 
-T < t < T, whereas when 1 < T < 2 we have T-2<r<2-T. 



B.l. Derivation of Eq. (4-. 4^) 



From Eq. (3.8) for the invariant measure we find, for two channels 

P ( - 2 \t 1 ,t 2 ) = C{tx-t 2 ) 2 . 
The normalization constant is found from 



1 = C 



(ti - r 2 ) 2 dridT 2 = — 



1 ,T 

dT dr- t 2 

ii J-T 



so that C=G. 

The distribution of the spinless conductance is thus 

rT' 



2 /-2-T 

dT dr- t 2 

J. JT-2 



c 

6' 



(B.2) 
(B.3) 



so that 



w {2) {T < 1) = 2T 



/ dT'5{T-T') I dr-r 2 + dT'5{T-T') / dr. 

JO J —T' Jl JT'-2 

w {2 \l < T < 2) = 2(2 — T) 3 . 



(B.4) 



(B.5) 



These two equations can be condensed in the result of Eq. ( |4.42[ ) 
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B.2. Derivation of Eq. U-4\ ) 

From Eq. ( |3.8| ) for the invariant measure we find, for two channels 

\P1T2 

The normalization constant is found from 
f = C 

2 r 





Tl - T 2 


ch~i 0IT2 






c 










[ dr 


2 


Jo 


J-T 



VT 2 - T 2 



dT 



2-T 



dr- 



2H 



T-2 



VT 2 - T 2 



~3°' 



(B.6) 



(B.7) 



so that (7=3/4. 

The distribution of the spinless conductance is thus 



8 7 



u> (1) (T) = iH/ dT'8(T-T') 



T' 



dr- 



2M 



,2 ,2-T' 

/ dT'J(T - T ) / 

Jl JT'-2 



2\t\ 



so that 



W W(T< 1) = ~T, 



;«(! < T< 2) = ~ 



r- 2vr^T 



which are Eqs. (4.41) we wanted to prove. 



(B. 



(B.9) 



B.3. Derivation of Eq. (4-44) 

The integration region selected by the 8 function in Eq. ( 4.39| ) is N ~ 1 dimensional. In the particular 
case T < 1, each of the t q 's (a = = 1, . . . , N) varies, over that region, in the interval < r a < T. This 
is clearly illustrated in Fig. Bl for TV = 2, where the relevant region is the segment extending from 
( T i: T 2) = (0, T) to (ti,t 2 ) = (T, 0). For T < 1 the corners of the TV-dimensional hypercube are not 
relevant; they become relevant for T > 1. We can thus write, for T < 1 



«#>(T<1) = C 



_ ^t(/3) 

'TV 



- 70 - 70 V a=l / a<b 



10) 



Introducing the new variables o~ a — T a /T (a — 1, • ■ -TV) which vary in the interval (0, 1), we have 



< ) (T<1) = C 



N 



f 1 ...[ 1 s(T-T£ t * a ) n k - * b f t n ^ n n 

70 ^ V a=l / a<b c 



(B.ll) 



0=1 



a<6 



which behaves as the power of T indicated in Eq. (4.44). 
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